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Chemical-potential dependence of pi and sigma meson masses is analyzed at both real and imaginary chem- 
ical potentials, and /ii, by using the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model that pos- 
sesses both the extended Z3 symmetry and chiral symmetry. In the fii region, the meson masses have the 
Roberge-Weiss periodicity. The dependence of the meson masses becomes stronger as temperature increases. 
We argue that meson masses and physical quantities in the /xr region will be determined from lattice QCD data 
on meson masses in the /ii region by using the PNJL model, if the data are measured in the future. 
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I. INTRODUCTION 

Extensive studies for exploring the phase diagram of Quan- 
tum Chromodynamics (QCD) have been done at finite temper- 
ature (T) and chemical potential (/i). In the study of the quark- 
gluon system, lattice QCD (LQCD) is a powerful method if 
11 = 0. LQCD, however, has the well known sign problem 
when the real part of fi is finite; for example, see Ref. [1] and 
references therein. Therefore, several approaches such as the 
re weighting method [,2J, the Taylor expansion method [3J, the 
analytic continuation from imaginary chemical potential /zi to 
real one /iR P, 's', ^ and so on are suggested to circumvent 
the difficulty, but those are still far from perfection. 

Constructing the effective model is an approach comple- 
mentary to the first-principle LQCD simulation. For example, 
the phase structure and light meson masses at finite T and 
/iR are extensively investigated by the Nambu-Jona-Lasinio 
(NJL) model S i, [M 111 and the the Polyakov-Ioop ex- 



oop ex- 

lijii 



tended Nambu-Jona-Lasinio (PNJL) model 1 12 ,113 , 

[IlSlillMllSlinillllllMlisilael 12711^ 

The NJL model can treat the chiral symmetry breaking, but 
does not possess the confinement mechanism. Meanwhile, the 
PNJL model is designed 1 14] to treat the confinement mecha- 
nism approximately in addition to the symmetry breaking. In 
this sense, the PNJL model is superior to the NJL model. 

In the PNJL model with two flavor quarks, the model pa- 
rameters are usually determined from the pion mass and the 
pion decay constant at T = /i = and LQCD data at finite T 
and zero /i; see Sec. II for the details. However, the strength of 
the vector-type four-point interaction can not be determined at 
zero ^ and then remains as a free parameter, although the loca- 
tion of the critical endpoint of the chiral phase transition in the 
/iR region is found to be sensitive to the strength 1 1 1, 28, 31]]. 
Thus, it is highly nontrivial how large the strength of the 
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vector-type interaction is and whether the PNJL model well 
simulates the jj, dependence of the phase structure and light 
meson masses. This should be tested directly from QCD. For- 
tunately, this is possible in the /ij region, since LQCD is feasi- 
ble there because of the absence of the sign problem. Further- 
more, it is possible to determine the strength of the interaction 
,for example the vector-type four quark and the higher-order 
multi-quark interaction, from LQCD data in the /ii region, as 
proposed by our previous paper [32]. 

In addition, the canonical partition function Zc{n) with 
real quark number n can be obtained as the Fourier transform 
of the grand-canonical one Zgc(^) with /i = i^i = i9T [33.] : 



(1) 



Thus, the thermodynamic potential, ^2qcd(^^) = 
-Thi(ZGc(^')), at finite 6 contains the QCD dynamics 
at real n and hence at finite /<r in principle. Therefore, we 
can confirm the reliability of the PNJL model in the /iR 
region by comparing the model results with LQCD ones in 
the /ii region. 

Roberge and Weiss (RW) found that fiQci){S) has a pe- 
riodicity, >!7qcd(^?) = ^qcb{0 + 27rfc/3), in the /ij region 
llssll . where k is any integer. The RW periodicity indicates 
that QCD is invariant under the extended Z3 transformation, 
that is, the combination of the Z3 transformation and the co- 
ordinate transformation 9 + 27r/3, as shown later in 
Sec. ini see our previous works (3^ for the details. At the 
present stage the PNJL model is only a realistic model that 
possesses both the extended Z3 symmetry and chiral symme- 
try. The PNJL model results are then consistent with LQCD 
ones particularly in the (^-dependence of the Polyakov loop, 
the quark number density and the chiral condensate 1132(1 . In 
the PNJL model, we do not need any extrapolation from the 
III to the /iR region, since the model calculation is feasible 
at finite /iR with the input parameters determined so as to re- 
produce LQCD data in the /xj region. We call this procedure 
the imaginary chemical-potential matching approach (the 9- 
matching approach) in this paper. 

In this paper, using the PNJL model, we predict the /i de- 
pendence of pi and sigma meson masses in the real and imag- 
inary /i regions, and argue that meson masses and physical 
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quantities in the /ir region will be determined from LQCD 
data on meson masses in the fii region by using the PNJL 
model, if the data are measured in the future. Concretely, 
the following four points are argued. First, we show that the 
meson masses have the RW periodicity. Second, the 9 de- 
pendence of meson masses is found to become large as T 
increases. We then recommend that LQCD simulations on 
meson masses in the /ij region be made at higher T near 
the pseudo-critical temperature Tc of the deconfinement phase 
transition at /i = 0. Third, the validity of two extrapolations 
from /i = /ii to /i = /ipi is discussed by comparing results of 
the extrapolations with the PNJL one. In recent LQCD cal- 
culations at finite 9, only small lattice sizes are taken, so that 
pion mass evaluated in Ref. ll34ll is unnaturally large. This 
indicates that the bare quark mass (mo) taken there is rather 
large. Finally, we discuss how sensitive meson masses are 
to the value of mp, comparing two cases of toq — 5.5 and 
80 MeV. This sort of model prediction is important before do- 
ing heavy LQCD calculations with large lattice size in the /ij 
region. 

We briefly explain the PNJL model in section II and present 
equations for meson masses in section III. In section IV, nu- 
merical results are shown on the /i-dependence of pi and 
sigma meson masses, and the validity of two extrapolation 
methods is discussed. Section V is devoted to summary. 



II. PNJL MODEL 

In this section, we briefly review the PNJL model; see |132 
for the details. The Lagrangian density of the two-flavor PNJL 
model is 

-CpNJL =q{i"fyD'' - mo)q + Gs[{qq)'^ + [qil^Tq)"^] 

~U{<P[AmAlT), (2) 

where = 9" + iA", q denotes the quark field with two 
flavor and the current quark mass toq. The field A" is defined 
as A" = (5o3^a%- with the gauge field A"^, the Gell-Mann 
matrix Aq and the gauge coupling g. The matrix r stands for 
the isospin matrix and Gs denotes the coupling constant of 
the scalar-type four-quark interaction. For simplicity, we ne- 
glect the vector-type four-quark interaction, since it does not 
change the conclusion of this paper qualitatively. 

In the PNJL model, the gauge field A^ is treated as a homo- 
geneous and static background field Aq, that is, A^ — Sq^Aq. 
In the Polyakov gauge, the Polyakov loop and its Hermitian 
conjugate, ^ and 3, are diagonal in color space: 



1 



7trc(i), 



with 



L = e'^-'^ = diag( 



(3) 



(4) 



where (fja and 0;, are classical variables and = lA^). 

We make the mean field approximation (MFA) to the quark- 
quark interactions in dU, as follows. In dU, the operator prod- 
uct qq is first divided into qq — a + (qq)' with the mean field 



a = (qq) and the fluctuation (qq)'. Ignoring the higher-order 
terms of (qq)' in the rewritten Lagrangian and resubstituting 
(qq)' = qq — a into the approximated Lagrangian, one can 
obtain a linearized Lagrangian based on MFA, 



MFA 
PNJL 



--qiij^d" + ijoAi - M)q 

-U{$[A13[A],T), 



(5) 



where M is the effective quark mass defined by M = toq — 
2Gs(J- In (|5]l, use has been made of {qijsfq) = 0, because 
the ground state is assumed to be invariant under the parity 
transformation. In the MFA Lagrangian 'CpNjl, quark fields 
interact only with the homogeneous and static back ground 
fields Aq and a. Hence, we can easily make the path integral 
over the quark field to get the thermodynamic potential per 
unit volume. 
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PNJL 



= -2Nf 



3 En 



(27r)3 
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P 

+ Um+U, 



-3/3(_Ep-^0l 



(6) 



where ^ = /iR + i^i = /iR + iT9, Ep = + AP and 
Um = Gsa\ 

The thermodynamic potential J7pnjl is invariant under the 
extended Z3 transformation 1I32I1 . 



■^{9)e'^. 



<P{9)t 



27rfc 
' 3 
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This is easily understood by introducing the modified 
Polyakov loop, •P' = e*^^ and 'F = e^*^^, invariant under the 
the extended Z3 transformation since i7pNjL is described 
as a function of only the extended Z3 invariant quantities, 'F, 

^, a and e^'" 



PNJL 



= - 2Nf 



d3p 



(27r)3 
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The physical quantities X = cr, and ^ are determined by 
the stationary conditions dQ'p^:ii^/ dX — 0. These equa- 
tions include the 6'-dependence only through the factor e'^'^, 
indicating that the X have the RW periodicity, X{9) = 
X{9 + 2nk/2>). Inserting the solutions back to i7pNjL, one 
can see that i7pNjL also has the RW periodicity, J?pnjl {9) = 
npN.iL{6 + 27rA:/3); see fs^l for the details. 

We take the three-dimensional momentum cutoff because 
this model is nonrenormahzable. 



(27r)3 
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(9) 
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Hence, the present model has three parameters mo, A and Gg. 
We take A = 0.6315 GeV and Gs = 5.498 GeV^^ so as to 
reproduce the pion decay constant f-^ — 93.3 MeV and the 
pion mass M-^ = 138 MeV at T = ji = 0, when a reahstic 
quark mass mg = 5.5 MeV is taken ifsillssll . 

We use U of Re f. OTtIi fi tted to LQCD data in the pure gauge 
theory at finite T BOW^ : 



U 



b2{T) 
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b2{T) = ao + ai{'^ 
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(10) 
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Following Ref. |32], we take Tb 190 MeV so as to repro- 
duce the pseudo-critical temperature Tc of the deconfinement 
phase transition at /x = evaluated by LQCD; specifically, Tc 
is 170 MeV in the PNJL model, whereas it is 173 ± 8 MeV 
in full LQCD calculations llssll . Thus, the present four param- 
eters are determined from the pion mass and the pion decay 
constant at T = /i = and LQCD data at T > and /i = 0. 
However, if the vector-type four-quark interaction (97^(7)^ is 
added to C, the strength can not be determined at /i = 0, since 
its mean field n = {q'joQ) is zero there. The strength will be 
determined by physical quantities at finite 6 such as meson 
masses, if they become available in the future. 



III. MESON MASS 

In this section, we consider pion and sigma meson and de- 
rive equations for the meson masses, following Ref |l2?l. Cor- 
relators of current operators carry physical mesons with the 
quantum number. The pseudoscalar isovector current with the 
same quantum number as pion is 



Jp^ix) = q{x)ij5T''q{x) 



(12) 



and the scalar isoscalar current with the same quantum num- 
ber as sigma meson is 



Jsix) = qix)q{x) ~ {q{x)q{x)). 



(13) 



The Fourier transform of the mesonic correlation function 
(0|r (^J^{x)JI\0)^ |0) is defined as 



C'lliq^) = i I d^a;e"?-"(0|T (jj (x) jf (0)) |0) 



2\xab 



(14) 



where £, = P for pion and S for sigma meson and T stands 
for the time-ordered product. Using the random phase ap- 
proximation (the ring approximation), one can obtain the 
Schwinger-Dyson equation for G^j at T = /i = where 

<P = <P = 0, 

= n^di^) + 2Gsn^di^)c^i (15) 

with the one-loop polarization function 



(-0 



(2^ 



TTir^iSip + q)r^lSiq)), (16) 



where S{q) is the quark propagator in the Hartree approxima- 
tion and — Fp = i^^r^ for pion and = Fs = 1 for 
sigma meson. In the random phase approximation, the solu- 
tion to ( fTSl ) is given as 



Gee = 



l-2Gsi7ee((?2)- 



(17) 



Noting that meson mass A/^ (^ = tt and a) is a pole mass of 
C^^{q^) and taking the rest frame q = (qo, 0), one can get an 
equation for as 



[l-2G«7I«(go)]|,„=M, =0- 
The explicit forms of 7Tpp(go) and 7Jss('?o) are given as 
^fpp(<7o) 



(18) 



= -iXr 

= 4i trc 
> 

nssiqo) 

= — iXr 
= 4i tr^ 



[ij5TiS{p + qo/2)ij5fiS{p - (7o/2)] 
d^P {p + qo/2){p-qo/2)-A'P 



(27r)4 
74 



(2^)4 [{p + qo/2r-AP+ie] 
1 



[{p ~ (?o/2)2 - AP + ie 



(19) 



0^[^S{p + qo/2)^S{p-qo/2)] 

d^p ip + qo/2)ip-qo/2) + AP 
(2^)4 [{p + qo/2r-AP+ze] 
1 



[{p - go/2)2 + M2 + ie 



(20) 



where Tr denotes troirac 5>) trs 



pin * 



) trc (X" trpia 



When T and 11 are finite, the corresponding equations are 
obtained by the replacement 



Po 



lUJn 
'4„ 



_(£p_ 
(2^ 



^ + = (2n + \)itT + 11 + iAi, 
d^p 



(27r)3 



(21) 



Also in this case, an equation for meson mass Af^ is of the 
same form as ( fTsT i, but the polarization functions are obtained 
in more comphcated forms: 

nMl^) - -2iVf [2^(m) - qlB{qo,y)] (22) 

for pion and 

i7ss(go) = -2Ni[2A{ii) - {ql - ^hP)B[qo,y)], (23) 

for si gma meson. Here, A{^) and B{qo, fi) denote loop inte- 
grals II39I1 defined by 



Ai^,) = 



B[qo,p) 



d^P - fpNJhjEp - fi) - fpNJhjEp + n) 
(27r)3 2Ep 

(24) 

d^p 1 - /pnjl(-Ep - /i) - /pnjl(£^p + m) 
(2^)3 Ep{ql^4Ep) 

(25) 
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where /pnjl {Epin) is the modified Fermi-Dirac distribution 
function [2 1.1 defined by 

/pNJl(£'p — fi) 

^ 1 -(_ 3(|(^ _|_ ^g-P{Ep-ti)^g-P{Ep-fi) _|_ g-3/3(£;p-Ai) ■ 

(26) 

The distribution function /pnjl (^'p + m) is obtained by re- 
placing -> ^ — > ^ and ^ ^ In actual PNJL cal- 
culations, the imaginary part of meson mass is assumed to 
be negligible in (fTSl i: this assumption is exactly satisfied when 
Mx is smaller than twice the dynamical quark mass, 2M, and 
approximately satisfied for slightly above 2M. 

In the case of fi — ifii = iTO, the corresponding distribu- 
tion functions 

hm^{Ep - ire) 

depend only on the extended Z3 invariant quantities, 'F, fp, a 
and e^*^. The distribution function /pnjl(£'p — iT9) is ob- 
tained by replacing 9 ~9, ^ W and ^ ^P. Therefore, 
the distribution functions are also extended Z3 invariant, and 
hence they have the RW periodicity. Furthermore, the polar- 
ization functions 77^^ depend on 9 only through the modified 
Fermi-Dirac distribution function /pnjl(-Ep ± iT9). There- 
fore, the meson masses have the RW periodicity. 

IV. NUMERICAL RESULTS 

First, we investigate the /i dependence of pi and sigma me- 
son masses, A/^ and Ma-, in the /iR and /ij regions. It is found 
from (l24l i and jZST l that A/^ and are symmetric under the 
interchange /i ^ — /i, indicating that they are functions of /i^. 
Figure [T] presents the fj? dependence of and A/^r in the 
case of T = 160 MeV that is near the pseudo-critical temper- 
ature Tc — 170 MeV of the deconfinement phase transition 
at /i = 0. They are smooth at = 0, as expected. This 
makes it possible the analytic continuation of meson masses, 
Af7r(/i) and Ma{^), from jj. = ifii to /i = /ip. This is im- 
portant for LQCD simulation, although the PNJL calculation 
does not need the analytic continuation. In the right-half panel 
representing the /ip region, M-^ and Ma agree with each other 
when ^0.1 GeV^. This clearly exhibits that the restoration 
of chiral symmetry takes place at >0.1 GeV^ ll2lll . In the 
left-half panel representing the /ij region, Afjr and M^ look 
almost constant in this scale, but one can see oscillations of 
the masses on closer inspection shown by the inset and Fig. |2] 
that presents M^^ and M^ as a function of 9. 

In our previous paper ||32|]. we showed that 9-odd quantities 
with the RW periodicity such as df2pj^ji^/d9 and the imagi- 
nary part of tZ/ are discontinuous at 9 = tt/3 (mod 2tt/3) and 
T > 190 MeV. This first-order transition is called the RW 
phase transition. Meanwhile, 6'-even quantities with the RW 
periodicity such as a and the real part of 'F are discontinuous 



in their derivative (3^. This indicates that the meson masses 
M^ have the same property as the chiral condensate. Actually, 
Fig.|2]shows that the A/j are 6'-even functions with the RW pe- 
riodicity and then not smooth on the RW phase transition line. 




-6.1 0.1 

[GeV^] 



Fig. 1: The ^'^ -dependence of sigma and pi meson masses, Ma and 
Mtt, at r = 160 MeV. The inset represents the pion mass near fi = 
0. 

As an important property, the phase of oscillation is op- 
posite to each other between M^ and M^^. This can be un- 
derstood from their slopes at fi^ — 0; dMa/dfi^ < while 
dM-^/diJ? > as shown in the inset of Fig. [T] These slopes 
reflect the chiral symmetry restoration appearing at /i^ > 0. 
Another point to be noted is the difference in their amplitudes; 
the amplitude of oscillation is relatively larger in M^ whereas 
very small in M-„. This can be understood from ( |22] | and ( l23T i; 
the latter has a term depending on the effective quark mass 
M and consequently reflects the chiral symmetry restoration 
directly, while the former does not. 

The meson mass M^{9) is a 6'-even function with the RW 
periodicity in the /ij region. This means that M^{9) can be 
expanded in terms of cos(3fc0) with integer k: 

M^{9) = ^afe(r)cos(3fc6l). (28) 

When T — 160 MeV, the normalized coefficients 
ah{T) / ao{T) are about 0.1 % for k — 1 and negligibly small 
for k > 2. The neglect of the normalized coefficients with 
fc > 2 is a good approximation at T near and below Tc. Par- 
ticularly in the strong coupling limit of LQCD, the ak with 
k > 2 are known to be zero li4Q1 . Hence, M^{9) is approxi- 
mated into 

Mf = A{T) cos(36i) + C(T), (29) 

where A{T) and C{T) are determined from M^ at 6* and 

7r/3 as 

A{T) ^ [M^{T, 9 = 0)- M^{T, 9 = 7r/3)]/2, (30) 

C(T) = [M^{T,9 = 0) + M^{T,9 = 7r/3)]/2. (31) 
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QI{T^li) 0.16 0.2 0.24 
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Fig. 2: 9 dependence of (a) sigma and (b) pi meson masses. Dotted 
curves denote the results of T = 160 MeV and solid curves corre- 
spond to tlie results of T = 200 MeV. Note that the scales of the 
vertical axes are different between panels (a) and (b). 



Figure [5] shows A and C for the case of A/jr. They 
smoothly increase with increase in T except a dip around 
T = 237 MeV that comes from a threshold effect due to 
TT — > quark + antiquark. The ratio A/C also becomes large 
as T increases. This means that the ^-dependence of M^^ be- 
comes stronger as T increases. It is then preferable that LQCD 
simulations will be made at higher temperature near Tc except 
the dip temperature, in order to determine the parameters of 
the PNJL model more accurately. The ^-dependence of me- 
son masses at lower T can be extracted from LQCD data at 
higher T by using the PNJL model. 

If A = 0, Ajfj will have no dependence, as shown by 
(|29] l. The fact that A/C is small, then, indicates that the 9 
dependence of is rather weak. Figure |4] presents the T- 
dependence of pion mass at 6* = and 7r/6, where results of 
9 = n/Q (9 = 0) are represented by solid (dotted) curves. 
As expected, the overall behavior is almost the same between 



Fig. 3: T-dependence of A(T) and C(T) for the case of pion mass. 



the two cases and the difference between the two is not large, 
although the crossing point where A/jr — 2Af is shifted to 
higher T as 9 increases. 

The extrapolation of from fj? < to /i^ > was neces- 
sary for LQCD at least so far Actually, LQCD data at /i^ < 
were extrapolated to the /i^ > region by assuming some 
fitting functions 14,5", '6']. In the PNJL model, physical quanti- 
ties are calculable directly in the /i^ > region, if the param- 
eter set is determined from LQCD data in the fi^ < region. 
Here we assume the present PNJL result as a result of this 
(^-matching approach and compare it with results of usual ex- 
trapolations in order to test the validity of the extrapolations. 
We consider two extrapolations. The first one is 

JVdn/T) = A{T) cosh(3Ai/T) + C{T). (32) 

When fi = i/ii, Eq. ^ equals to Eq. This form of M^' 
has an exponential /i-dependence at /i = fi^. This extrapola- 
tion is then called the exponential extrapolation in this paper; 
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Fig. 4: T-dependence of sigma and pion masses at ^ = and n/6. 
Twice the constituent quark mass M is also shown. 



a similar form is used in Ref. JH] for thermodynamical quanti- 
ties rather than meson masses. The second is the polynomial 
extrapolation [4^^ in which a polynomial of /i^ is taken up to 
eighth order 

In Fig.|6] Af^ calculated with the exponential and the poly- 
nomial extrapolation are compared with the result of the PNJL 
model. In the two simple extrapolations, their parameters are 
fitted to the PNJL result for /^^ < q at T = 160 MeV. At 
/i^ < 0, M^* agrees with Mr, calculated with the PNJL model 
within thickness of curves. This means that the Ofe for k > 2 
are tiny in the expansion (l28T l. The two simple extrapolations 
give almost the same result for /i^ > 0. Thus, the polyno- 
mial extrapolation often used so far is essentially equal to the 
exponential one. The PNJL result and the results of the two 
extrapolations coincide accurately up to /xr/T ~ 1. 

Next, we discuss mathematically on the the analytic contin- 
uation of A/j from /x — ifii to = /iR + i/ij. First, we assume 
that LQCD gives Af^ numerically in the region. In order to 
make the analytic continuation, we need an analytic form of 
that is valid in the /ij region. The Fourier expansion se- 
ries (l28T l gives such an exact form of Af^, since A/^ is a 6'-even 
function with the RW periodicity. As shown below, A/^ has no 
9 dependence in the low-T limit and then all the except ao 
tend to zero in the limit. The partition function Zqc {S) with 
a finite value of 9 is equivalent to Zqq (0) with the boundary 
condition q{x, 1/^) = — exp {i9)q{x, 0) for the quark field 
q lHH. In the low-T limit where a period 1/T of the imagi- 
nary time becomes infinite, the value of Zqc{0) does not de- 
pend on how to take the boundary condition and then has no 
9 dependence. The PNJL model can reproduce this property, 
as proven in Ref. js^l . 

As a result of the analytic continuation from fi — ifii to 

= Mr + the series becomes 



Fig. 5: /^^ -dependence of M-^ and Af^* at T=160 MeV. The solid 
curve represents the PNJL result, the dotted curve does a fitted one 
by a polynomial up to (fi^)^ and the dashed curve does a fitted one 



This is a natural extension of the exponential extrapolation 
( l32b . This form is valid, if the series converges. The condition 
for the convergence is 



R — lim 

k — >oo 



ak+i{T) cosh{3{k + l)n/T) 



ak(T) cosh(3fc/i/T) 



< 1. 



(34) 



In general, R depends on /i and T. Now we consider 
the case of fi — In the low-T limit, all the a/, ex- 

cept ao tend to zero, as mentioned above, while the fac- 
tor cosh (3(fc + 1)/^r/T)/ cosh (3fc/iR/T) diverges for all k. 
Thus, there is a possibility that R is less than 1 even in the low- 
T limit. However, it is impossible to evaluate the extremely 
small coefficients (k > 1) from LQCD data in the finite 9 
region. Also for large T near Tc, the PNJL calculation shows 
that the coefficients ak with k > 2 are negligibly small. It 
is then difficult to evaluate the small coefficients from LQCD 
data with finite errors in the finite 9 region. Hence we propose 
to use the ^-matching approach based on the PNJL model in- 
stead of the exponential and polynomial extrapolations and 
the analytic continuation ( [33] l. 

Finally, we check an influence of the current quark mass 
mo on the /i-dependence of A/jr, since large mo is taken in 
LQCD simulations at finite 9. The result adopting mo = 80 
MeV is shown in Fig. |6] for example. This figure indicates 
that, with heavy mo, the sign of dAI.^/d^j?' at ^j? = becomes 
opposite to the case of light mo shown in Fig. [T] This should 
be noticed. 



V. SUMMARY 



M^{^i) =Y^ak{T) cosh(3A:/i/T). 

A;=0 



(33) 



We have analyzed, by using the PNJL model, the /i- 
dependence of pi and sigma meson masses in both the /iR 
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Fig. 6: -dependence of pion mass at T = 160 MeV in the case of 
mo = 80 MeV. Definition of curves is the same as in Fig. ??. 



and Hi regions. In the /ij region, the meson masses M^{9) 
{£, — TT and a) are even functions of 9 with the RW periodic- 
ity: M^{e) = M^{-0) = M^{e + 2fc7r/3). This property is 
the same as that of the chiral condensate. The RW periodic- 
ity indicates that in general M^{6) are oscillating with 9. The 
amplitude of the oscillation becomes large as T increases. We 
then recommend that LQCD calculations be done at higher T 
near the pseudo-critical temperature Tc of the deconfinement 
phase transition at /i — , in order to determine the param- 
eters of the PNJL model more accurately. The 6'-dependence 
of meson masses at lower T can be extracted from the LQCD 



data at higher T by using the PNJL model. As for pion mass 
Mt^{9), it should be noticed that the phase of the oscillation is 
rather sensitive to the value of the current quark mass mg. 

It is possible to do the PNJL calculation with a parame- 
ter set common between the /iR and /ii regions. Hence, we 
can argue that meson masses and other physical quantities in 
the /xr region can be extracted from LQCD data on meson 
mass at (1) finite T and zero ji and (2) finite T and nonzero 
/ii, by using the PNJL calculation the parameters of which 
are fixed to the LQCD data. The present PNJL model has 
three parameters, too. A, Gg, in the quark sector and one 
parameter Tq in the gauge sector. As an extension of this 
minimal PNJL model, the vector-type four-quark interaction 
and/or the scalar-type eight-quark one are occasionally added 
to the present Lagrangian i32[l4lll . All the parameters can be 
determined from LQCD data at the two cases; a trial is shown 
in Ref. |41]. The polynomial and exponential extrapolations 
used so far are accurate at /ir/T < 1. The 0-matching ap- 
proach based on the PNJL model is expected to be one of the 
most reliable methods to get physical quantities at finite /iR, 
if LQCD data on M^{6) become available in the near future. 
We recommend that LQCD simulations be done systemati- 
cally with the same lattice size between two cases of (1) finite 
T and zero /i and (2) finite T and nonzero /ij. 
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